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Abstract: This paper summarizes liiglilights of the OGl session of the 30th International Cosmic Ray 
Conference, held in Merida (Yucatan, Mexico). The subsessions (OGl.l, OGI.2, OGI.3, OG1.4 and 
OGl. 5) summarized here were mainly devoted to direct measurements, acceleration and propagation of 
cosmic rays. 
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The OGl session was splitted in 5 subsessions, 
devoted to direct measurements of cosmic rays 
by balloons and satellites (OGl.l), cosmic ray 
composition (OGl. 2), cosmic ray propagation 
(OGl. 3), cosmic ray acceleration (OG1.4) and in- 
trumentations and new projects (OG1.5). The 
number of papers discussed in the OGl session 
gives a feeling of the huge number of results pre- 
sented during the conference: there were 138 pa- 
pers presented (66 of which were oral presenta- 
tions). About 37 of the oral presentations reported 
on observational results, while about 29 reported 
on theoretical or phenomenological work. It is ob- 
vious that a meaningful summary of this work im- 
plies a selection of highlights. An apology is due 
to all those that will not see their work properly 
discussed here. It should also be understood that 
in order to avoid doing a simple list of the results 
it is needed to put things in context, and this often 
requires weighing the results according to the opin- 
ions of the writer, which do not coincide necessar- 
ily with the opinions of the Community at large. I 
hope that this may be a starting point for further 
discussion on the topics that I wiU touch upon be- 
low. 



Introduction 

The origin of cosmic rays is a problem that has 
been haunting scientists for almost one century 
now. Solving this scientific problem means putting 
together numerous pieces of a complex puzzle, in 
which the acceleration processes, the inner dynam- 
ics of the sources, the propagation, the chemical 
composition all fit together to provide a satisfac- 
tory and self-consistent global picture. We are not 
there yet, though an increasingly larger number of 
pieces are finding their place in the puzzle. The In- 
ternational Cosmic Ray Conference is a privileged 
place to observe the puzzle taking shape. 

This rapporteur paper is organized as follows: in 
§3 I present a subjective bird's eye view of how 
things stand in the field, in order to make it easier 
for the reader to put in the right context the ex- 
tremely wide range of results presented at the Con- 
ference. In §4 I summarize some important obser- 
vational results about spectra and chemical com- 
position of cosmic rays (including the electronic 
component) that have been presented. In §5 I il- 
lustrate some recent developments in the under- 
standing of the acceleration of cosmic rays. In 
§6 I discuss the problem of relating the observed 
cosmic rays to the sources through propagation in 
the Galaxy. The transition from the galactic com- 
ponent to the extragalactic cosmic ray component 
is briefly discussed in §7 together with some is- 
sues related to the propagation of ultra high energy 
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cosmic rays in the Galaxy and in the intergalactic 
medium. I conclude in §8. 

A bird's eye view on cosmic rays 

The best known property of cosmic rays is their all- 
particle spectrum. Most experiments agree, at least 
qualitatively, that the spectrum consists of four re- 
gions: 1) at low energies (below ~ 10 GeV) the 
observed spectrum is flat as it is affected by so- 
lar modulation. 2) For energies lOGcV < E < 
3 X 10 ^^eV, the spectrum can be fit with a power 
law with slope 2.7. 3) At energies between 
3 X lO^^GcV and ^ lO^^cV the slope grows to 
~ 3.1. 4) The upper boundary of region 3 is not 
well defined, so I will take some liberties in defin- 
ing region 4. In terms of Physics this reflects our 
uncertainties in defining the region of transition to 
extragalactic cosmic rays (see §7). 

The chemical composition of cosmic rays is a cru- 
cial piece of information and is becoming weU de- 
termined at energies below the knee, where direct 
measurements can still be carried out, while it is 
more uncertain and model dependent at higher en- 
ergies. The composition of low energy cosmic rays 
provides important hints to the acceleration pro- 
cesses and the propagation of cosmic rays through 
the interstellar medium (ISM). Especially impor- 
tant in this respect are the abundances and spectra 
of elements such as Boron, Beryllium and Lithium, 
which are mainly produced as secondaries of pri- 
mary cosmic rays. The ratio of secondary to pri- 
mary (for instance B/C) cosmic ray fluxes provides 
a unique tool to characterize the diffusion proper- 
ties of the ISM. Existing measurements of this ra- 
tio as a function of energy suggest that the diffu- 
sion coefficient scales with energy as D{E) oc E", 
with a w 0.6, at least at rigidities below ~ 10 GV, 
while it is not clear whether at higher energies the 
slope remains constant or there is a flattening. 

Elements such as Ge and Ga can potentially allow 
us to discriminate between volatility based and first 
ionization potential based acceleration processes, 
and therefore provide information about the domi- 
nant acceleration sites. Other elements (e.g. ^^Co 
and ^^Ni) also provide us with precious informa- 
tion on the mean time between the production of 



the material which is accelerated and the actual 
time when it is accelerated to cosmic ray energies. 

Current direct measurements are filling the gap 
to the knee region, thereby providing a way to 
match and cross-check the indirect measurements 
of spectrum and composition carried out by ob- 
serving and modelling extensive air showers at en- 
ergies across and above the knee. The data col- 
lected by the KASCADE experiment suggest that 
the knee in the all-particle spectrum may be an 
artifact of the superposition of sharper knees in 
the single components, mainly in the light chem- 
icals, such as H and He. The proton spectrum 
measured by KASCADE has a pronounced knee 
at energy ^ 10^ GeV. A knee in the He compo- 
nent is also seen, but the fluxes become more un- 
certain while moving to larger masses. It is not 
known as yet what could be the explanation for the 
knees in these components, though it is plausible 
that they may reflect inefficiency of the accelerator. 
In a strictly rigidity-dependent approach, a knee in 
the iron spectrum could be expected at ^ 3 x 10^ 
GeV. This simple argument opens the way to an 
equally simple but far reaching impUcation: the 
spectrum of Galactic cosmic rays should end at en- 
ergies around ~ 10^'' eV, so that cosmic rays of 
larger energies must be accelerated in extragalactic 
sources. Unfortunately the KASCADE results are 
not fully consistent with those of some other exper- 
iments, especially the Tibet array. The main diffi- 
culty of all ground based experiments, including 
those operating in the ultra high energy range, is 
the poor understanding of cosmic ray interactions 
in the atmosphere, which may be very problematic 
to infer basic information as the chemical compo- 
sition of the primary particles. 

Three sets of measurements will help us figure 
out whether the physical picture suggested by the 
KASCADE data is correct (assuming that the data 
themselves find proper confirmation): 1) direct 
measurements of the cosmic ray spectra and chem- 
ical composition should extend as far as possible 
towards the knee and possibly across it; 2) addi- 
tional data and rehable models for the description 
of the shower developments (possibly checked ver- 
sus future LHC data) are crucial at energies around 
and above the knee; 3) the chemical composition 
in the energy region between lO^'^ eV and 10^^ eV 
needs to be measured reliably. 
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Fortunately, important steps ahead in all these di- 
rections are being done and some important results 
will be discussed below. 

What about the sources? The paradigm based on 
supernova remnants (SNRs) as the main sources 
of galactic cosmic rays remains the most plau- 
sible, but the smoking gun that would turn the 
paradigm into a well established theory is still 
missing. The most important news in this direc- 
tion is represented by the recent observations of 
SNRs in gamma rays by Cherenkov imaging tele- 
scopes (see [31] and references therein for a re- 
cent review) and the high resolution observations 
of the X-ray emission from the rims of several 
SNRs (see [59] and references therein for a recent 
review). The former have presented us with some 
evidences of TeV gamma ray emission, possibly of 
hadronic origin. The latter have provided us with 
strong evidence of efficient magnetic field amplifi- 
cation, which is in turn required in order to reach 
the high energies observed in cosmic rays, and are 
a consequence of efficient cosmic ray acceleration. 
From the theoretical point of view, there have been 
several new developments in our understanding of 
the mechanism of diffusive shock acceleration in 
SNRs but also related to the propagation of cos- 
mic rays both in the Galaxy and in the intergalac- 
tic medium. 1 will discuss some of these issues in 
more detail below. 

New IMeasurements 

The direct measurements of the flux of cosmic 
rays and of their chemical composition, carried 
out by using balloons and satellites, has always 
played a crucial role in advancing our understand- 
ing of both acceleration and propagation of cos- 
mic rays, at least at energies below the knee. At 
higher energies the low fluxes due to the steeply 
falling spectrum make it necessary to use ground 
arrays which observe cosmic ray induced air show- 
ers. The two techniques are clearly complementary 
and one of the biggest problems has always been 
to cross-calibrate the two. For the first time this 
goal seems to be at least in sight, in that the di- 
rect measurements are extending to ~ 10 — 100 
TeV, therefore approaching the knee region. Nu- 
merous new results on spectra of different chem- 



ical elements have been presented, from balloon 
flights (CREAM, ATIC, Tracer, TIGER, Bess- 
Polar, PPB-BETS), from satelhtes (preliminary re- 
sults from PAMELA) and even some shower ex- 
periments (Tibet Array and HESS). The overlap 
between all these techniques in the energy region 
around the knee appears to be of crucial impor- 
tance, especially to unveil the origin of the knee 
in the cosmic ray spectrum. 

Balloons and Satellites 

The CREAM Collaboration has presented impres- 
sive results collected during a total of 70 days of 
flight of their balloon experiment. The data were 
collected during the record breaking CREAM-1 
flight (42 days between 12/16/05 and 01/27/05) 
and a second flight (CREAM-11) lasted 28 days 
(between 12/16/05 and 01/13/06). Convincing evi- 
dence was presented of the excellent charge resolu- 
tion of the experiment {Aq ^ 0.2 electron charge), 
precious tool to reconstruct the spectra of different 
chemicals (see [50, 23] for a technical discussion). 

In Fig. 1 I reproduce the plots as shown by [55] 
(see also [63]) illustrating the spectra of protons 
and heUum nuclei as measured by CREAM-I (red 
circles). The CREAM spectra are compared with 
those of other experiments as hsted in the figure 
with different symbols. The spectrum of He nuclei 
as measured by CREAM is compatible with that 
measured by ATlC-2, but appears to be somewhat 
flatter than that measured by ATIC-1. 

The spectra of Carbon and Oxygen measured by 
CREAM-II were presented by [65] and [2]. The 
results on the C/O ratio confirm the primary na- 
ture of both nuclear species [51]. An overview 
of the CREAM results and future developments (a 
third and fourth flight) were presented by [55] . The 
spectra of the chemicals presented by the CREAM 
collaboration do not show appreciable differences 
in the slope, with the possible exception of the he- 
lium spectrum that might be slightly harder than 
the proton spectrum. In fact this trend is probably 
present even in the ATIC-2 spectrum. In Fig. 1 
the shaded area is supposed to show the range of 
fluxes as measured by ground arrays. The width 
of the shaded region provides an estimate of the 
uncertainties due to either the interaction models 
adopted for the analysis or systematics in the dif- 
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Figure 1 : Spectra of protons (left panel) and helium nuclei (right panel) as measured by CREAM-I (figure 
from [55]) compared with the results of other experiments. 



ferent experiments. An important point that this 
figure shows is that finally direct and indirect mea- 
surements are starting to overlap in the knee re- 
gion. This will be of great importance to under- 
stand the reason for the appearance of the knee in 
the all-particle cosmic ray spectrum. 

The ATIC balloon experiment has presented im- 
pressive results on the B/C ratio and mass compo- 
sition of cosmic rays below the knee. In fact the 
ratios N/0 and C/O were also measured. The B/C 
ratio, as I discuss in §6, is a crucial indicator of the 
mode of cosmic ray propagation in the interstel- 
lar medium. Low energy measurements (E < 30 
GeV) carried out by HEAO-3 [30] show that the 
B/C ratio scales with rigidity as ^ usually 
interpreted as a proof of the primary (secondary) 
origin of Carbon (Boron) nuclei. 

The ratio as measured by ATIC-2 is reported in Fig. 

2 [49] together with the previous results of HEAO- 

3 and the predictions of two theoretical models 
(lines) (see [49] for some discussion and refer- 
ences). Despite the apparent flattening of the mea- 
sured B/C ratio in the ATIC data, the error bars are 
still too large to infer solid conclusions on this is- 
sue. 

The all-particle spectrum and the average mass 
number as a function of energy as measured by 
ATIC-1 are shown in Fig. 3 [1] compared with the 
results of other experiments, as indicated in the fig- 
ure. There is a substantial agreement among all the 
data sets shown. At least in part the small offsets 
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Figure 2: B/C ratio as a function of energy per nu- 
cleon as measured by ATIC-2 (squares) [49] and 
compared with the results of HEAO-3 and the pre- 
dictions of two theoretical models (lines). 
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may be explained in terms of systemaic errors in 
the energy determination, amplified by the multi- 
plication by E^-^. 

The mean logarithmic mass number found by 
ATIC in the energy range 10^ — lO^GeV shows a 
general trend to a heavier composition and appears 
to match well with the results of RUNJOB, CASA- 
BLANCA, DICE and KASCADE in the knee re- 
gion. At energies above the knee (not reached by 
ATIC) the other experiments show quite different 
trends. This prevents us from reaching a satisfac- 
tory explanation of the origin of the knee and of 
the transition from galactic to extragalactic cosmic 
rays (see §7). 

Very accurate measurements of the flux of different 
chemical elements below the knee were presented 
by the TRACER Collaboration [21] and are shown 
in Fig. 4 (filled dots) for nuclei between Oxy- 
gen and Iron. The fluxes are multiplied by differ- 
ent normalization factors to make the plot clearer, 
and they are compared with the results of HEAO-3 
and CRN. No appreciable difference between the 
slopes of the spectra of these nuclei was detected, 
all slopes being around ~ 2.7. 

The TRACER Collaboration also presented the 
B/C ratio in the energy region around ^ 100 
GeV/amu, being fully consistent with other mea- 
surements. The predicted error bars on the B/C 
ratio with a 30 days balloon flight of TRACER 
should allow to finally pin down the slope of the 
ratio as a function of energy up to ~ 1 TeV/amu. 
While indicators such as the B/C ratio provide 
us with important information on the propagation 
of cosmic rays in the Galaxy, the abundances of 
super-heavy elements (Z > 30) tell us about the 
acceleration regions, though their flux is more than 
~ 1000 times smaller than the Iron flux. The 
flux of these elements has been measured during 
two balloon flights by the TIGER Collaboration in 
2001 and 2003 (a total of about 50 days of observa- 
tions). The fluxes in the energy region between 0.3 
and 10 GeV/nucleon are reported in Fig. 5 (from 
[54]). The lines represent the predictions of two 
models for acceleration, one based on first ioniza- 
tion potential (FIP) and the second on volatility, but 
both with a standard solar system composition of 
the ambient medium. The abundances of ^^Ga and 
^^Ge appear to be inconsistent (if taken at the same 
time) with both theoretical models. One should 
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Figure 4: Differential spectra of chemical elements 
from Oxygen to Iron as measured by Tracer [21] 
and compared with previous results from CRN and 
HEAO-3. 



however recall that a standard solar system com- 
position near the accelerator is all but granted. 

Important information on the origin of cosmic rays 
and their propagation can also be indirectly gath- 
ered by observing electrons. Their propagation in 
the Galaxy is affected by diffusion, as for ions, 
and energy losses (mainly synchrotron and inverse 
Compton scattering losses). At energies larger than 
20-100 GeV the loss time due to ICS and syn- 
chrotron becomes shorter than the diffusive escape 
time from the Galaxy, so that the spectrum steep- 
ens by one power compared with the source spec- 
trum. At energies larger than ~ 1 TeV the flux 
of electrons at Earth is expected (and observed) to 
decline, presumably as a result of the discreteness 
in the spatial distribution of the sources. In this 
sense several groups have been engaged in a search 
for a possible signal of the closest sources of elec- 
trons around the solar system, the most promising 
sources being Vela, the cygnus loop and Mono- 
gem. The expected signal would consist of a bump 
in the diffuse flux of electrons at energies in excess 
of ^ 1 TeV. The proximity of the source would 
also result in a small anisotropy in the arrival di- 
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Figure 3: All particle spectrum (left panel) and mean logarithmic mass number (right panel) as measured 
by ATIC [1]. 
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Figure 5: Relative abundances of heavy elements 
as observed by the TIGER experiment [54]. 



Figure 6: Spectrum of electrons (multiplied by E^) 
as observed by PPB-BETS [64] compared with the 
results of previous measurements. 

rections. The spectrum and anisotropics have been 
presented by PPB-BETS. 

To date, the spectrum [64] (Fig. 6) does not show 
evidence for the possible appearance of a nearby 
electron source. Even the anisotropy is claimed to 
be fully consistent with an isotropic distribution of 
arrival directions [64] . 

The BESS-Polar experiment has presented several 
interesting results, ranging from low energy cos- 
mic ray spectra to spectra of antiprotons and lim- 
its to the flux of anti-helium nuclei. The flux of 
cosmic ray protons in the energy region below 
^ 10 — 30 GeV is heavily affected by solar mod- 
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Figure 7: Spectram of protons as observed by 
BESS- 1997 (empty circles) and BESS-Polar (filled 
circles) [32]. The difference between the two spec- 
tra is due to the different levels of solar modulation 
at the time of the measurements. 

ulation. In Fig. 7 1 report the results presented by 
[32]: the thick dots are the recent results of BESS- 
Polar, compared with those of BESS- 1997 (open 
circles), at the time of the solar minimum. It is ev- 
ident the effect of the enhanced solar modulation 
while the sun gradually approaches the next so- 
lar minimum. The upper solid line represents the 
spectrum of protons in the interstellar medium as 
inferred from BESS-98 measurements. 

The propagation of cosmic ray nuclei in the in- 
terstellar medium also results in the production 
of many secondary products (nuclei, antiprotons, 
electron-positron pairs, gamma rays, neutrinos). 
Accurate measurements of these secondaries, es- 
pecially antiprotons and positrons is instrumental 
to the search for faint signals, such as those coming 
from the annihilation of non-baryonic dark mat- 
ter. The flux of anti-protons and the ratio of anti- 
protons to protons in the energy region below few 
GeV as measured by BESS-Polar is reported in 
Fig. 8. In the left panel, the antiproton spectrum is 
compared with some theoretical predictions of the 
same quantities for some models (solid line: propa- 
gation with GALPROP; dotted line: standard leaky 
box with solar modulation; dash-dotted line: black 
hole evaporation). In the right panel, the ratio an- 



tiproton/proton flux is compared with predictions 
of a drift model with different tilt angles of the so- 
lar magnetic field. The thick dots and empty circles 
refer again to BESS-Polar and BESS- 1997 respec- 
tively. In the right panel the effect of a different so- 
lar modulation in the two periods is evident again. 
As expected, this effect almost completely disap- 
pears in the ratio of the antiproton/proton flux. The 
fluxes of antiprotons observed by BESS -1997 and 
BESS-Polar are both fully consistent with the one 
expected on the basis of a propagation model with 
solar modulation, thereby implying significant up- 
per limits to the flux of antiprotons from, for in- 
stance, neutralino annihilation in the Galactic dark 
matter halo. BESS has also focused on the detec- 
tion of anti-nuclei. The quest for why our universe 
is (almost) completely made of matter instead of a 
mixture of matter and anti-matter is a fundamen- 
tal one, as Nature is expected to have produced a 
(almost) symmetrical universe. Though it is likely 
that the small excess of matter over anti-matter is 
what now makes most of the present universe, the 
search for islands of antimatter or traces of anti- 
nuclei has never quite finished. The limit imposed 
by BESS-Polar and the limit that will possibly be 
reached by the next BESS flight are shown in Fig. 
9, together with the results of previous measure- 
ments. 

The measurement of the cosmic ray spectrum and 
chemical composition up to energies of 1 TeV, 
the search for dark matter, anti-matter and ex- 
otic particles, and finally bits of solar and mag- 
netospheric physics are among the goals of the 
PAMELA satellite, which has been succesfully 
launched on June 15, 2006. The instruments on- 
board are collecting data, but only preliminary re- 
sults were presented at the conference. A review 
of PAMELA preliminary results and expectations 
has been presented in [52]. Particularly interest- 
ing, and suggestive of the quality of the results that 
we should expect, is the spectrum of protons and 
helium, reported in Fig. 10. This preliminary anal- 
ysis does not show evidence for any difference in 
the slopes of the spectra of H and He in the energy 
region below ~ 500 GeV. In both cases the best fit 
power law has a slope of 2.73. 

Several simulated results were also presented by 
the PAMELA Collaboration, mainly aimed at 
showing the potential for discovery of indirect sig- 
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Figure 8: Left Panel: Antiproton spectram measured in BESS-Polar compared with BESS(95+97) results. 
The lines illustrate predictions of different theoretical scenarios. The lines that pass through the data points 
are the predictions of standard propagation calculations (GALPROP and leaky box models). Right Panel: 
p/p ratio measured by BESS(95+97) and BESS-Polar I (2004) compared with drift model calculations at 
various Solai^ magnetic field tilt angles. 




Figure 10: Preliminary spectra of protons and helium from PAMELA. 
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Figure 9: Limits on the flux of anti-Helium ob- 
tained by BESS-Polar [32]. 

nals of dark matter annihilation (in the form of 
anomalous features in the antiproton and positron 
spectra), and the limits achievable in the search for 
anti-matter. The limit on the ratio He/ He is ex- 
pected to improve by about 1 order of magnitude 
in about three years of PAMELA operation with 
respect to that achieved by BESS-Polar (see Fig. 
9). 

PAMELA is expected to play a crucial role in the 
identification of spectral features in the positron 
and antiproton spectra as they can possibly result 
from dark matter annihilation, provided the mass 
of the dark matter candidate particle is in the right 
energy range (below 500 — 1000 GeV). 

Ground Experiments 

Since most results from ground experiments were 
presented in other sessions of the ICRC, the talks 
on this subject in the OGl session were sparse. 

The HESS Collaboration has presented the re- 
sults of a very interesting detection of the direct 
Cherenkov (DC) light from cosmic rays impacting 
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Figure 1 1 : Spectrum of iron nuclei (in fact nuclei 
with Z > 24 obtained by HESS using the detection 
of DC light. The HESS results (with S YBILL and 
QGSJET as models for interactions) are compared 
with the results of other measurements. 

the atmosphere. This technique, proposed as a pos- 
sible tool to measure the chemical composition of 
cosmic rays above 10 TeV by [58], has been used 
by HESS to measure the flux of heavy nuclei (Fe- 
like). The direct Cherenkov light is produced by 
the primary nucleus while entering the atmosphere 
and before the first interaction that generates the 
shower High up in the atmosphere the density 
of air is lower and therefore the Cherenkov cone 
is more collimated with respect to the Cherenkov 
cone of the light generated by lower energy parti- 
cles in the shower. As a result the total Cherenkov 
emission from a shower is expected to have a broad 
region with a hot pixel corresponding to the DC 
light. This emission has in fact been success- 
fully detected by HESS. Since the intensity of the 
Cherenkov signal scales with the square of the 
charge of the parent nucleus, it is best to seach for 
the signal from Fe nuclei provided the energy does 
not exceed a few hundred TeV, in order to avoid 
that the shower emission overshines the DC light. 
The spectrum of Iron nuclei as measured by HESS 
is reported in Fig. 11 [22] where it is compared 
with the results of other experiments. The two sets 
of data points refer to QGSJET and SIBYLL as in- 
teraction models. 

The search for the DC light is also being pursued 
by VERITAS [62], while dedicated instruments are 
being planned (e.g. TrICE [61]). 
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The HESS Collaboration has also presented the 
preUminary results of their measurement of the 
electron spectrum [29], consistent with measure- 
ments carried out by other experiments within the 
systematic uncertainties. 

A relevant contribution to the understanding of the 
nature of the knee in the all-particle spectrum has 
been provided by Tibet Array. The Collaboration 
has presented the most recent measurements of the 
proton and helium spectra in the knee region, and 
the fraction of heavy nuclei as a function of energy, 
in the range 10^ < E < lO'^ GeV. 
The measurements of the chemical abundances in 
this energy region are still controversial. The Tibet 
data on the proton spectrum appear in rough agree- 
ment with those of RUNJOB and JACEE (both 
with QGSJET and SIBYLL as interaction mod- 
els) but in rather apparent contradiction with the 
KASCADE data, especially when the KASCADE 
data are analyzed using SIBYLL. The latter show a 
steeper spectrum of protons and a somewhat higher 
normalization of the flux. The discrepancies are 
even more evident in the case of the helium spec- 
trum. The Tibet results on heUum roughly match 
with those of RUNJOB, but not with data from 
JACEE and KASCADE. The absolute flux of he- 
Uum as measured by KASCADE is almost one or- 
der of magnitude higher than that measured by Ti- 
bet Array. The situation in the knee region contin- 
ues to be confused and this prevents the possibility 
of vmveiUng the origin of the knee. 

Acceleration of Cosmic Rays 

The main mechanism for the acceleration of cos- 
mic rays (not only galactic) remains diffusive par- 
ticle acceleration at shocks. Many contributions 
presented at the conference have focused on the 
determination of the efficiency of particle acceler- 
ation, the dynamical effect of the accelerated parti- 
cles and the magnetic field ampUfication generated 
by the accelerated particles through streaming in- 
stability or due to turbulent amplification. 

I will start by discussing some issues related to 
the acceleration of Galactic cosmic rays in SNRs, 
though some of the conclusions will be of wider 
appUcabiUty. 



One of the most recent advancements in the theory 
of particle acceleration at non-relativistic shocks, 
especially for SNR shocks, has been the calcula- 
tion of the dynamical reaction of the accelerated 
particles onto the background plasma. The effect 
has been discussed already in the 80's in the con- 
text of two-fluid models (see [44] for a review), 
and later addressed numerically by direct solution 
of the time dependent equation of diffusive tran- 
port coupled with the equations of mass, momen- 
tum and energy conservation of the background 
plasma (see for instance [11]). In the late 90' s a sta- 
tionary semi-analytical solution of this set of equa- 
tions was found in the form of an integral equation 
[43] for a spatially constant diffusion coefficient. A 
general semi-analytical method was recently pro- 
posed by [8], vaUd for any choice of the diffusion 
coefficient. 

The structure of the shock is changed by the accel- 
erated particles due to the pressure they exert on 
the background plasma, though in a collisionless 
manner. For ordinary diffusion coefficients, which 
are increasing functions of particle momentum, 
higher energy particles can travel further away than 
low energy particles: a fluid element approach- 
ing the shock surface therefore feels an increasing 
pressure due to accelerated particles and as a con- 
sequence it slows down. This leads to the forma- 
tion of a precursor upstream of the shock (which is 
now called subshock) in which the fluid velocity, as 
seen in the (sub)shock reference frame decreases 
while approaching the shock. The effective com- 
pression factor felt by particles crossing the shock 
and precursor is now a function of momentum, be- 
ing bigger at large momenta and smaller at low 
momenta. The immediate consequence is that the 
spectrum of particles accelerated at shocks is not a 
power law, as expected in test-particle theory, but 
rather a concave function of momentum, harder at 
high momenta and softer at low momenta. This 
nonlinear system is found to self-regulate itself and 
to be able to reach relatively high efficiency of par- 
ticle acceleration. 

This complex nonlinear system has however ad- 
ditional interesting aspects, concerning the maxi- 
mum momentum that the particles can be accel- 
erated to. As was recognized first in [39, 40], 
the maximum momentum achieveable in SNRs is 
exceedingly low compared with the knee energy 
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unless the diffusion coefficient is Bohm-like and 
magnetic field amplification takes place. Even in 
this case the maximum energy falls short of the 
knee energy by a factor ^ 100, unless the ampli- 
fication turns strongly nonlinear, namely 5B/B ^ 
100. An investigation of the process of particle ac- 
celeration at cosmic ray modified shocks with self- 
generation of strong turbulence was discussed in 
[9, 18] and presented at this Conference by [20]. 

The most important clue in this field has however 
come in the last few years from high resolution X- 
ray observations of the rims of several SNRs. The 
thickness of the rims is related to the loss length 
of high energy electrons radiating X-rays by syn- 
chrotron emission. The observed brightness pro- 
files lead to estimates of the downstream magnetic 
field of the order of a few hundreds iiG [60]. These 
estimates were confirmed and presented here by 
[16, 15] together with their implications for mul- 
tifrequency observations of specific SNRs. 

These calculations are based on a numerical time- 
dependent solution of the problem of particle ac- 
celeration at modified shocks. The strength of the 
magnetic field is not self-consistently calculated 
but rather fit to the observations at a given time 
(the observation time) and not evolved in time. The 
calculations of [15] show that the multifrequency 
data for supernova RXJ17 13.7-3946 are explained 
rather well from radio to X-rays and gamma rays 
(Fig. 12 from [15]). In particular the gamma ray 
data from HESS are interpreted as the result of 
pion production and decay. The magnetic field re- 
quired for the fit is 126/i(j, again consistent with 
the general trend of strong magnetic field ampli- 
fication observed in other rermiants. The strong 
magnetic field causes the electron density to be 
low, thereby reducing the flux of ICS gamma rays. 
Similar consequences can be inferred for the RX 
J0852.0-4622 (Vela Jr.) [16]. It is also worth 
noticing, in Fig. 12 (and others that have been 
presented for other SNRs) that both the predicted 
radio-X and gamma ray spectra show a clear indi- 
cation of a curvature (concavity). This characteris- 
tic, as discussed above, is one of the most distinc- 
tive consequences of particle acceleration at mod- 
ified shocks. Signatures like this should be care- 
fully looked for in future, better data at all wave- 
lengths. In particular, future gamma ray observa- 
tions with GLAST might allow us to confirm or 



reject the hadronic origin of the gamma ray emis- 
sion. 

Two models have been discussed for the amplified 
magnetic field observed in several SNRs. Quanti- 
tatively they potentially lead to the same level of 
magnetization, but the principle is different in the 
two cases. 

The traditional explanation is based on stream- 
ing instability induced by the accelerated particles. 
The instability leads to resonant growth of Alfven 
modes [10] and typically SB/B ~ 20 - 30 [9]. 
By using the results of quasi-linear theory for the 
growth rates, [20] showed how to determine the 
diffusion coefficient and implement it in a self- 
consistent solution of the problem of acceleration 
at modified shocks. At the beginning of the Sedov 
phase this type of instabihty leads to reach approx- 
imately proton energies of 10^ GeV, while at 
later stages of the SNR evolution, the maximum 
momentum decreases in a way that depends on 
how the magnetic field amplification drops. 

In [12] it was proposed that non resonant modes 
may grow faster than resonant modes in some cir- 
cumstances. The analysis of Bell was based on a 
MHD treatment of the background plasma and the 
non resonant modes are often considered as strictly 
related to this assumption. The presentation of [19] 
has demonstrated that the same dispersion relation 
is obtained in the context of a purely kinetic ap- 
proach, if particle acceleration is efficient, as it 
is expected at cosmic ray modified shocks. The 
imaginary and real parts of the frequency of the 
propagating waves are plotted in Fig. 13 in units 
of /crLfi as a function of the wavenumber k in 
units of 1/ri 0- Here o is the Larmor radius of 
the particles in the background magnetic field and 
Vs is the shock velocity, taken here as lO^cm s~^. 
The efficiency in acceleration was assumed to be 
r] = UcR/pVg =0.1. The peak of the growth rate 
if not at krL,o ^ 1 which is again a confirmation 
that the faster growing modes are non resonant. 
However the same authors of [19] also showed that 
the peak moves towards fcr^.o ~ 1 when either the 
shock velocity decreases or ij decreases. This sug- 
gests that at later epoques during the SNR evolu- 
tion the resonant mode may become dominant. 

The search for the non resonant modes found by 
[12] was also carried out by using Particle-in-Cell 
(PIC) simulations. The results of such studies were 
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Figure 13: Real and Imaginary part of the fre- 
quency of the waves excited by the streaming of 
cosmic rays at a shock with velocity lO^cm s^^ 
and efficiency of particle acceleration 77 = 0.1. For 
these parameters the peak of the growth rate oc- 
curs at fcr^ 3> 1, showing the non resonant nature 
of the instability [19]. 



presented by [47]: no evidence was found of the 
fastly growing Bell mode in their simulations. It 
remains to be seen whether this result is due to the 
anomalous values of the parameters adopted by the 
authors (for instance the large density of non ther- 
mal particles needed to satisfy the condition of effi- 
cient acceleration and the artificial value of the ra- 
tio rrie/mp), or if actually shows that the instability 
may be suppressed due to some physical mecha- 
nism yet to be identified. 

It is important to stress that streaming instabil- 
ity generates the magnetic field amplification up- 
stream of the shock front. The turbulent field is 
then advected downstream and the perpendicular 
components are compressed at the shock surface. 

A second model for the generation of strong field 
was discussed by [37]: the model is based on 
the (realistic) assumption that density perturba- 
tions Sp/p ^ 1 are present upstream (see Fig. 14 
for a schematic view). The density perturbations 
induce corrugations in the shock structure which 
in turn produce turbulent eddies that twist the field 
lines frozen in the plasma leading to magnetic field 
amplification. The numerical simulations that the 
authors illustrated show that SB/B ^ Ma, where 
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Warped Shock Front 




Figure 14: Schematic view of the turbulent mag- 
netic field amplification scenario proposed by [37]. 

Ma is the Alfven Mach number. In conditions 
which are typical of SNRs one may easily ex- 
pect SB/B ^ 100 - 1000 downsti-eam of the 
shock. The effects of such amplified magnetic field 
on the acceleration of particles need some further 
discussion: if the field in the upstream region is 
mainly parallel to the shock normal, and there is no 
streaming instability (as assumed in [37]) then the 
acceleration process is inefficient, being weakly af- 
fected by the downstream strong fields alone. In 
fact the acceleration time is the sum of the diffu- 
sion times upstream and downstream. If the field is 
only amplified downstream, the acceleration time 
is dominated by the upstream region, and it re- 
mains too long to lead to efficient acceleration. 

The problem can be avoided by assuming that the 
shock is perpendicular In this case the upstream 
acceleration time is substantially shorter because 
the perpendicular diffusion coefficient is smaller 
than the parallel one [36], thereby mimicking am- 
plification of the magnetic field from the point of 
view of acceleration. For perpendicular shocks the 
amplification by streaming instability is expected 
to be suppressed, therefore the two mechanisms 
could be considered as complementary and might 
take place together though with different efficien- 
cies as related to the relative orientation of the 
shock normal and the ambient magnetic field. 

Efficient particle acceleration at SNR shocks mod- 
ified by the dynamical reaction of cosmic rays, to- 
gether with the amplification of the magnetic field 
provide a consistent picture of the origin of cos- 
mic rays: protons in the SNR environment can be 



accelerated to ^ 10^ GeV, while elements with 
higher charge (Z) may be accelerated to energies 
which are Z times larger. The maximum ener- 
gies are reached at the beginning of the Sedov 
phase while at later stages the maximum energy 
decreases, mainly because the magnetic field am- 
plification is less efficient (at least for amplification 
due to streaming instability). The spectrum of cos- 
mic rays observed at the Earth is the superposition 
of the cosmic rays trapped behind the shock and 
released after adiabatic decompression and the par- 
ticles which may leave the system from upstream 
due to lack of confinement in the accelerator 

The calculations presented by [14] represent the 
first attempt to convolve the results from particle 
acceleration with given diffusion properties of the 
interstellar medium to determine the all-particle 
spectrum of Galactic cosmic rays. Good agree- 
ment with data is obtained after summing over all 
chemical elements, though there are numerous pa- 
rameters that can be tuned to obtain the fit. The 
most important point is however that the end of 
the Galactic spectrum is expected to be at energies 
around a few 10^^ eV. I will discuss this point in 
more detail in §7. 

Our understanding of cosmic ray acceleration at 
SNR shocks has certainly improved considerably 
in the last few years, as also shown by the im- 
pressive quality of the fits to multifrequency obser- 
vations of SNRs from radio to gamma ray wave- 
lengths. However, it may be useful to briefly dis- 
cuss some aspects of the problem that are not so 
well understood and that are nevertheless crucial to 
explain observations and to have a consistent pic- 
ture of the origin of cosmic rays. I will limit myself 
to mentioning three of such aspects. 

• 1) A general trend of the non linear theories 
of particles acceleration at modified shocks 
is to lead to very large total compression 
factors, incompatible with the values (~ 7) 
which provide a good fit to data. It is usu- 
ally believed that this reduction may be at- 
tributed to turbulent heating in the precursor 
[45, 13]. However to date we have no reli- 
able theory of how to treat this phenomenon, 
initially put forward to avoid magnetic field 
amplification to SB/B ^ 1 and currently 
used also in cases where the amplification is 
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strongly non linear. One can adopt recipes 
that qualitatively work, but do not ensure the 
correctness of the quantitative results. 

• 2) The flux of particles that may escape the 
shock region from upstream of a modified 
shock can be an appreciable fraction of the 
energy flux. We do not know as yet how to 
carry out a detailed calculation of this flux, 
despite the fact that cosmic rays detected at 
the Earth might be dominated by this com- 
ponent rather than by the ones trapped be- 
hind the shock. The spectrum of particles 
that escape a SNR at each time during the 
Sedov phase is expected to be very peaked 
at energies around the maximum momentum 
achievable at that specific time. As proposed 
by [53] the observed flux might be due to an 
overlap of these peaked spectra integrated on 
the SNR evolution. 

• 3) As discussed above, the crucial element 
which makes cosmic ray acceleration pos- 
sible in SNR shocks is the magnetic field 
amplification. If this process does not take 
place, SNRs can hardly be the sources of 
Galactic cosmic rays. The amplification may 
however take place in different ways (res- 
onant and non resonant streaming instabil- 
ity, firehose instability, turbulent amplifica- 
tion, and possibly others). Each amplifica- 
tion process leads to a different saturation 
level and different scaling laws with the age 
of the SNR. More specifically, at different 
times in the SNR evolution, different ampli- 
fication processes may operate. This makes 
the calculation of the actual spectrum of CRs 
from SNRs highly non trivial. 

It is remarkable that all the work done and pre- 
sented at the Conference on shock modifica- 
tion and its phenomenological consequences was 
mainly concentrated upon non relativistic shocks 
and more specifically SNR shocks. 

A paper was presented [48] in which the authors 
stress the difficulty of accelerating particles at rela- 
tivistic shocks (in the test particle approximation), 
which reflects into very steep spectra. This find- 
ing adds to the fact, which is becoming increas- 
ingly more recognized, that the dependence on am- 
bient conditions is very important for the case of 



relativistic shocks. For instance, even the simple 
compression of the perpendicular component of a 
turbulent magnetic field crossing a shock surface 
may induce a substantial steepening of the spec- 
trum in the relativistic case (see also [42]), while 
being unimportant for non relativistic shocks. 

Propagation of Cosmic Rays 

The propagation of cosmic rays in the Galaxy is 
dominated by diffusion (possibly anisotropic) and 
advection if there is a wind blowing outwards. 
QuaUtatively it is very easy to illustrate what might 
be expected. Let us assume that the sources in- 
ject cosmic rays into the ISM with a rate Qp{E) oc 
E~^, where the index p stays for primary. At the 
leaky box level, the equation which describes the 
stationary density of cosmic rays in the Galaxy is 



Qp{E) = 



ripjE) 
r{E) ' 



where t{E) is the time needed for escaping the 
Galaxy. If the escape is solely due to diffusion and 
the diffusion coefficient scales as D{E) oc i?", 
then t{E) oc E~°'. It immediately follows that 
np{E) oc Qp{E)t{E) oc ^-(«+t). During propa- 
gation in the ISM secondary nuclei are produced at 
a rate Qs{E) « np{E)cmHv{p) oc £;-("+'>') (as- 
suming v{p) c). It follows, again at the leaky 
box level of approximation, that the equilibrium 
density of secondaries is ns{E) « Qs{E)t{E) oc 
^-(a+27) xhej-gfore the ratio of the secondary to 
primary equilibrium densities is Us/up oc E~°' oc 
l/D{E). This is the case of the B/C ratio dis- 
cussed in §4: measuring the energy dependence 
of the B/C ratio we can infer the diffusion coeffi- 
cient, or more in general the escape time as a func- 
tion of energy, which scales as 1 /D{E) if diffusion 
is the only process responsible for escape while it 
is somewhat more complex if there are other pro- 
cesses involved (for instance advection in a wind). 

As discussed in §4 the most recent measurements 
of the B/C ratio by CREAM, ATIC and TRACER 
extended to energies of 100 — 1000 GeV/nucleon, 
thereby providing us with a unique opportunity 
to understand what is the rate of escape of cos- 
mic rays as a function of energy right below the 
knee. It is worth recaUing that the data points 
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with small error bars at energies smaller than 30 
GeV/nucleon (see Fig. 2) return a slope of the B/C 
ratio ^0.6. The higher energy data presented by 
the three Collaborations have too large error bars 
so far to imderstand if such a slope remains un- 
changed. On the other hand one should recall that 
if indeed cosmic rays escaped from the Galaxy pro- 
portionally to E~°-^, and normalizing the escaping 
time in the 10 GeV/nucleon region, one would in- 
fer a too large anisotropy of cosmic rays around the 
knee [33]. It is therefore Ukely that somewhere be- 
low the knee the behaviour of the escape time with 
energy changes to a flatter behaviour. On the other 
hand, if this transition exists then one should ex- 
pect the appearance of a feature in the all-particle 
spectrum, which does not seem to be there. The 
problem remains open. 

The importance played by diffusion in the inves- 
tigation of the origin of cosmic rays is also clear 
from the wealth of presentations at the Conference, 
concerning different ways to treat diffusion. The 
standard way to describe diffusion of cosmic rays 
in the Galaxy has become the GALPROP code, 
which was used also in the work presented by 
[46] in order to calculate the elemental abundan- 
cies throughout the periodic table. GALPROP is 
also used to determine spectrum and spatial distri- 
bution of the secondary emissions, especially radio 
and gamma, which provide information on the dis- 
tribution of magnetic field and gas respectively. 

Some random-walk-like approaches to diffusive 
motion have been presented in [27, 34] and then 
applied for specific purposes such as the calcula- 
tion of the diffuse galactic gamma ray background 
or the gamma ray emission from diffuse sources 
such as the HESS source in the galactic bulge. 

The simple arguments reproduced above, which 
illustrate the basic aspects of the leaky box ap- 
proaches, are a wild oversimplification of a phe- 
nomenon which is in fact very complex. Below I 
will discuss some of these complications and how 
to describe at least some of them. 

In general, diffusion acts on top of a regular mo- 
tion of charged particles moving in a large scale 
magnetic field and is due to a turbulent compo- 
nent SB. In the case of the Galaxy the large scale 
field has a complex structure, made of spiral arms 
and a poorly known halo. Diffusion parallel to 
the background magnetic field is faster than dif- 



fusion in the direction perpendicular to the ordered 
field. The difference between the diffusion coef- 
ficients in the two directions (parallel and perpen- 
dicular) decreases when the level of turbulent field 
increases. 

The parallel diffusion coefficient can be evalu- 
ated in the context of quasi-linear theory provided 
6B/B <C 1. In the non linear regime there is 
no definite theory for the determination of the dif- 
fusion coefficients and numerical simulations be- 
come invaluable tools. In the weakly non linear 
regime analyical techinques can still be applied and 
lead to rather unexpected results: for instance a 
similar theoretical approach presented by [56] and 
consisting of a weakly non linear approach to par- 
allel diffusion, results in a scaling of the parallel 
diffusion coefficient with energy as D\\ oc E^-^ 
despite the fact that the spectrum of turbulence is 
Kolmogorov-like. This might have an important 
impact on the interpretation of the observed slope 
in the B/C ratio, though this slope would still re- 
main incompatible with the observed anisotropy of 
cosmic rays at the knee. 

For typical magnetic fields of few /xG, as in the 
Galaxy, numerical simulations can be used to de- 
rive the propagation properties of cosmic rays 
down to energies of 10^* — 10^^ eV. The results 
of one such investigation were presented by [25] 
and revealed several interesting new aspects of the 
problem. The simulation consists of propagating 
a large sample of charged particles and determine 
their escape time from toy models of the magnetic 
field of the Galaxy. At the same time, in order to 
better understand and interpret the results the au- 
thors also calculate the diffusion coefficients (par- 
allel and perpendicular to the direction of the reg- 
ular magnetic field). 

The simulations cover the range of turbulence 
strength 0.5 < SB/B < 2 always assumed to 
be distributed according with a Kolmogorov spec- 
trum, oc k~^/^. The authors find fliat flie 
energy dependence of the perpendicular and paral- 
lel diffusion coefficients is different: D± oc E", 
with a w 0.5 - 0.6 and oc E'^/^ (basically 
the same slope obtained from quasi-linear theory 
despite the non linearity). 

The results of the calculations of the escape times 
from the Galaxy clearly illustrate the diffuculty 
of the problem. The authors investigated several 
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toy geometries of the regular magnetic field of the 
Galaxy, starting with a purely azimuthal field, from 
which particle escape can only take place through 
diffusion perpendicular to the field (and therefore 
to the disc). Even in this simple case, the escape 
time is affected not only by (perpendicular) dif- 
fusion, but also by drifts induced by gradients in 
the regular field. The gradients may be simply the 
ones associated with the curvature of field lines but 
could also be due to gradients along the radial di- 
rection in the disc, and along the z-axis perpendic- 
ular to the disc. The full set of results are discussed 
in [25, 26]. 

In Fig. 151 reproduce the escape times (top panel) 
and the grammage traversed by cosmic rays (bot- 
tom panel) as presented by [25] for the simple case 
of a purely azimuthal regular field. A few inter- 
esting results are apparent: first, the escape time 
obtained with all values of 5B/B have a slope 
~ 0.6, consistent with the fact that the escape time 
in this field configuration is oc /D±, where H 
is the thickness of the magnetized halo. Second, 
at -E ~ 10^^ eV the escape time is appreciably 
affected by the drift induced by the curvature of 
the regular magnetic field lines (solid line). Third, 
the grammage inferred from the simulation is of 
order 0.5 — 3g cm^^ at 10^^ eV. If extrapolated 
to 10 GeV with a slope 0.6 this would lead to 
exceedingly large grammage (or equivalently too 
long escape times) in the energy region where this 
parameter can be inferred from B/C ratio and is 
~ 20f; cm^^. A similar problem, though at higher 
energies, was previously found by [66]. Many 
other cases discussed by the authors confirm the 
generality of these few conclusions and show that 
the problem of propagation has stiU many aspects 
which are poorly understood and can hardly be 
included in phenomenological approaches includ- 
ing simple ones such as the leaky box models or 
more complex ones such as GALPROP. The prop- 
agation of ultra-high energy cosmic rays in the in- 
tergalactic space is in a way simpler to describe, 
though a major source of uncertainties is due to 
our ignorance of the magnetic field (strength and 
topology) possibly existing between the sources 
and our Galaxy. An appreciable intergalactic mag- 
netic field can affect both the spectrum and the 
anisotropy of cosmic rays. While the overall shape 
of the spectrum is left almost unchanged, because 
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Figure 15: Escape times (top panel) and grammage 
(bottom panel) for cosmic rays propagating in a toy 
model of the Galaxy [25]. 
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of the universality shown by [7], some features 
may appear in the observed spectrum because of 
the so-called magnetic horizon [6, 41]: at the low- 
est end of the cosmic ray spectrum, the propagation 
time from the nearest source may exceed the age of 
the universe thereby suppressing the flux of cosmic 
rays. For magnetic fields of less than ~ 1 nG, this 
effect typically appears below 10^^ eV (the details 
depend on the topology of the field and the diffu- 
sion properties of the particles). At higher energies 
the presence of the magnetic field mainly affects 
the anisotropy (especially on small scales) in the 
arrival directions. In fact at energies > 10^^^ eV 
the flux received at Earth is expected to come from 
nearby sources, and the arrival directions to point 
back to the sources. However there has been a lot 
of discussion on whether this is actually the case 
(see for instance [57] and [28] for two rather dif- 
ferent results), because of the uncertainties on the 
strength of the magnetic field. A paper was pre- 
sented at the Conference [38] in which the propa- 
gation in a structured magnetic field was discussed 
in the assumption that some fraction of the energy 
in the hydrodynamic flows leading to large scale 
structure formation is converted to magnetic fields. 
The authors reach a conclusion which is somewhat 
in between those of [57] and [28]: only some of 
the UHECRs should point back to their source. All 
these results should however be taken with a grain 
of salt, the uncertainties in the determination of the 
magnetic field being very large. From the observa- 
tional point of view, if a correlation was found be- 
tween arrival directions of UHECRs and a class of 
sources, this would be the best proof of the weak- 
ness of the magnetic fields in the intervening inter- 
galactic medium. It would be harder to conclude 
one way or the other in case of absence of correla- 
tions. 

The uncertainties in the galactic magnetic field are 
also rather disturbing and affect appreciably our 
ability to obtain realistic results on the propagation 
of cosmic rays. In [35] the authors proposed to 
use data on the polarization of the CMB photons 
to infer information on the large scale structure of 
the galactic magnetic field. Such a goal would cer- 
tainly justify the effort. 



Where do Cosmic Rays become extra- 
galactic? 

The theoretical arguments illustrated in §5 suggest 
that if cosmic rays are accelerated in SNR shocks 
then the maximum energy of nuclei of charge Z 
should be Z X 10^^ eV. This conclusion is based 
on the evidences for magnetic field amplification at 
SNR shocks, as inferred from X-ray observations. 
A similar conclusion could however be reached 
based on direct observation of the proton spectrum 
by the KASCADE experiment. The proton spec- 
trum as measured by the Tibet-Array experiment 
at energies above the knee is harder than that mea- 
sured by KASCADE, but the two agree on the fact 
that at ^ 10^^ eV there is a softening of the proton 
spectrum. Moreover, both agree, at least qualita- 
tively, on the fact that the chemical composition 
becomes heavier above the knee in the aU-particle 
spectrum. 

On the observational side this situation suggests 
that future accurate measurements of the chemical 
composition at energies above the knee will be cru- 
cial to solve the problem of the origin of cosmic 
rays. From the theoretical side, the most striking 
conclusion is that the galactic component of cos- 
mic rays should end at around ~ a few lO^'^ eV 
and in this energy region be iron dominated. This 
is in clear contradiction with the traditional ankle 
model of the transition, which postulates that the 
transition takes place at ^ 5 x 10^^ eV as a result of 
the intersection of a steep galactic spectrum and a 
flatter extragalactic (proton dominated) spectrum. 

Recently two other models have been proposed, 
both of which locate the transition region at ener- 
gies IQi'^ - 10^® eV, the dip model [17] and the 
mixed composition model [3]. They are described 
at length in the rapporteur paper on session HE and 
in the Review paper by V.S. Berezinsky, therefore 
here I will hmit myself to a short description. 
In the dip model [5] the extragalactic spectrum is 
proton dominated (not more than 15% of He nuclei 
are allowed). The effective injected spectrum at the 
extragalactic sources is required to be E^^i with 
7 = 2.6 — 2.7. The dip appears as a feature due 
to pair production losses of protons on the cosmic 
microwave background radiation and its location 
in energy is very well defined and model indepen- 
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dent. The shape of the dip is also very weakly de- 
pendent upon most parameters, with the exception 
of the chemical composition of the extragalactic 
component, as mentioned above. The extragalac- 
tic spectrum extends down to lower energies, and 
at the position of the so-called second knee it is 
predicted to flatten. This is the region of the tran- 
sition from a steep galactic to a flatter extragalac- 
tic cosmic ray component. The chemical compo- 
sition in the transition region is predicted to suffer 
a sharp transition from iron dominated to proton 
dominated [4], and the transition is completed at 
energy ~ 10^^ eV. 

In the mixed composition model the extragalactic 
component corresponds to a flatter injection spec- 
trum ^ E~'^-^, with a composition which is a mix- 
ture of different chemicals with abundances that 
are roughly comparable with the source composi- 
tions inferred for SNRs but with much freedom in 
this respect. The propagation of these components 
from the sources to Earth leads to the formation 
of a complex spectrum that may fit the all particle 
spectrum fairly well. The transition from galactic 
to extragalactic cosmic rays in this model is more 
gradual and is completed at energy ~ 3 x 10^^ 
eV. The chemical composition is expected to vary 
slowly in the transition region, shifting from iron 
dominated (at low energies) to graduaUy Ughter to- 
wards higher energies. 

The main discriminating factor between the two 
models is the chemical composition in the transi- 
tion region. Present measurements do not allow to 
make this discrimination as yet. 
In the two models illustrated above, the Galactic 
cosmic ray component is obtained by subtraction 
of the predicted extragalactic flux from the mea- 
sured all particle spectrum. In the presentation by 
[24] the authors have tried a different approach, 
namely that of calculating the galactic component 
by assuming a source spectrum and using GAL- 
PROP to descrive propagation in the ISM. The 
main goal was to show that in general, for both 
models, the superposition of the galactic and ex- 
tragalactic components leads to the appearance of 
strong features in the aU-particle spectrum. The au- 
thors claim that the features are more evident in the 
mixed composition scenario. This result should of 
course be taken as an interesting suggestion to be 
further investigated. We do not have enough in- 



formation about the shape of the spectrum of the 
different chemicals and their maximum energies in 
order to run GALPROP and obtain results which 
can be interpreted in a unambiguous way. 

Conclusions 

The collection of high quality data is the main drive 
for the field of cosmic ray research. In the very 
low energy region (below ~ 10 GeV) the accu- 
rate measurements of the cosmic ray flux by BESS 
has, among other things, allowed to achieve an ex- 
cellent understanding of the effects of solar mod- 
ulation. At energies below the knee in the all- 
particle spectrum, several experiments (CREAM, 
ATIC, Tracer) are providing us with large statis- 
tics of data and correspondingly excellent quality 
spectra of nuclei spanning the all periodic table. 
The domain of ultra-heavy nuclei is being covered 
by the TIGER experiment. These direct measure- 
ments are gradually approaching the energy range 
around the knee thereby providing us with a test of 
consistency (or inconsistency) of data collected in- 
directly by ground experiments using atmospheric 
showers induced by the cosmic ray interactions. 
One of the goals of the research in this field in the 
years to come should be (and in fact it is already) 
to match the results of direct and indirect measure- 
ments of cosmic ray fluxes and chemical composi- 
tion. Thinking of the possibility of ultra-long dura- 
tion balloon flights to reach closer to or across the 
knee appears as a worthy effort. 

The achievement of accurate measurements of cos- 
mic ray spectra is also opening the way to the 
search for weak signatures of rare phenomena, 
such as propagation of antimatter from nearby re- 
gions of the universe and annihilation of dark mat- 
ter in the Galaxy. The latter has profound im- 
plications on the spectra of positrons and antipro- 
tons especially, besides creating features in gamma 
rays and at other wavelengths. The launch of the 
PAMELA satelhte in 2006 has been a milestone in 
this direction, and hopefully next ICRC will have 
plenty of presentations of exciting new results. 

The picture that emerges from direct cosmic ray 
measurements can be summarized as follows. 
There is a satisfactory agreement on the all-particle 
spectrum as measured by different experiments. 
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The spectra of the chemical species do not show 
statistically significant differences in the slopes. 
The abimdances of Ga and Ge do not seem im- 
mediately compatible with acceleration scenarios 
based on pure first ionization potential or volatil- 
ity. Different determinations of the B/C ratio are 
consistent with each other up to the highest mea- 
sured energy (~ 1 TeV), but the error bars are still 
large enough to leave open the possibility of a flat- 
tening in the slope of the B/C ratio as a function 
of energy. A long duration flight, for instance of 
Tracer, is likely to settle this issue. We recall that 
from the theoretical point of view a naive extrapo- 
lation of the observed B/C ratio to the knee region 
would lead to exceedingly large anisotropy. 

A few results of measurements carried out with 
ground experiments, such as HESS and the Tibet 
hybrid detector were presented in the OG session. 
The HESS collaboration presented the positive de- 
tection of the direct Cherenkov radiation, and used 
it to infer the spectrum of iron nuclei (more cor- 
rectly of nuclei with charge Z > 24) with energy 
up to ~ 100 TeV. The measured spectrum appears 
to be in good agreement with other results in the 
same energy region. 

The Tibet Collaboration discussed the results of 
a measurement of the proton and heUum primary 
spectra as obtained with the Tibet hybrid experi- 
ment. The discrepancy between these results and 
those of the KASCADE experiment are clear, espe- 
cially for the helium spectrum. As stressed above, 
one way of facing these problems is to extend the 
direct measurements with high statistics of events 
to the knee region. 

The calibration of the ground arrays by an over- 
lap with direct measurements is a crucial goal to 
pursue, not only to understand the origin of the 
knee but also to describe correctly the transition 
from galactic to extragalactic cosmic rays. The two 
problems are tightly related to each other. 

There have been numerous presentations at the 
Conference on the theory of acceleration of cos- 
mic rays. From the phenomenological point of 
view it is clear that several independent pieces sup- 
port the possibility that the bulk of galactic cosmic 
rays are accelerated in SNRs: 1) X-ray observa- 
tions have shown that effective magnetic field am- 
plification at SNR shocks takes place; 2) multifre- 
quency observations of SNR RXJ1713 and Vela ju- 



nior are best fit if the observed gamma ray emission 
is of hadronic origin; 3) the spectra of radio and 
gamma rays show some evidence, though proba- 
bly not conclusive, of curvature (concavity), a phe- 
nomenon which is expected in non Unear theories 
of particle acceleration at shock waves, when the 
dynamical reaction of the particles is not negligi- 
ble, namely when acceleration is very efficient. 

From the theoretical point of view, the origin of 
the magnetic field amplification is being widely in- 
vestigated. Two mechanisms might be at work at 
the same time: streaming instability induced by 
the efficiently accelerated cosmic rays and turbu- 
lent amplification in a shocked medium with den- 
sity inhomogeneities. Other instabiUties, such as 
firehose instability, might also be at work. There 
have been several contributed papers on the growth 
rate of the non resonant streaming instability in 
the context of a kinetic approach, and as simulated 
with the help of PIC simulations. The latter do not 
seem to confirm the large growth rates predicted by 
quasi-linear theory, but additional work is needed 
to make sure that the initial conditions are set in a 
physically meaningful way. 

Efficient magnetic field amplification, when it is 
due to streaming instability, also imphes efficient 
particle acceleration, which in turn induces a con- 
cavity in the spectra of both the accelerated par- 
ticles and the radiations produced by them. Effi- 
cient magnetic field ampUfication is also required 
in order to explain cosmic ray energies at least in 
the knee region for protons and correspondingly 
higher for nuclei with higher charge. 

AU these ingredients are put together in multifre- 
quency investigations of SNRs, where evidence for 
magnetic field amplification comes from X-ray as- 
tronomy. Convincing evidence was presented that 
such multifrequency spectra, from radio to gamma 
rays, can be explained in the framework of parti- 
cle acceleration at modified shocks. It was also 
claimed that the propagated spectra of different 
chemical components, when summed together pro- 
vide a good fit to the observed all-particle spectrum 
of cosmic rays. 

Unfortunately propagation in the Galaxy is all but 
well understood, mainly because of our ignorance 
of both the regular and turbulent components of the 
galactic magnetic field: the competition between 
parallel and perpendicular diffusion, the presence 
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of a wind, the topology of the magnetic field in 
the spiral arms and in between the arms all affect 
the escape times and their dependence on energy. 
A hint of what might be going on could come from 
future extensions of the primary to secondary ratios 
to higher energies, possibly approaching the knee 
region. 

Acknowledgements 

I am grateful to the organizing Committee of 
the 30th Intemational Cosmic Ray Conference for 
their help. I am also grateful to the numerous sci- 
entists in Merida that 1 had the pleasure to talk to 
in order to better prepare this contribution. Finally 
I wish to express my gratitude to all my closer col- 
laborators, R. Aloisio, E. Amato, V. Berezinsky, D. 
Caprioh, D. De Marco, S. Gabici, G. Morhno, M. 
Vietri for continuous stimulating discussions. 

References 

[1] H. S. Ahnetal. All-Particle Spectrum Mea- 
sured by ATlC-1. In ICRC, Merida, Mexico, 
ICRC.OG.n73, 2007. 

[2J H. S. Ahn et al. Elemental Spectra from the 
CREAM-1 FHght. In ICRC, Merida, Mexico, 
ICRC.OG.1055, 2007. 

[3 J Parizot E. Olinto A. V. Kahn E. Allard, D. and 
S. Goriely. UHE nuclei propagation and the 
interpretation of the ankle in the cosmic-ray 
spectrum. Astron. and Astrop., 443:29, 2005. 

[4] Berezinsky V. Blasi P. Aloisio, R. and 
S. Ostapchenko. Signatures of the transition 
from galactic to extragalactic cosmic rays. 
Submitted to Phys. Rev. D (arXiv0706.2834), 
2007. 

[5 J Berezinsky V. Blasi P. Gazizov A. Grig- 
orieva S. Aloisio, R. and B. Hnatyk. A dip 
in the UHECR spectrum and the transition 
from galactic to extragalactic cosmic rays. 
Astropart. Phys., 27:76, 2007. 

[6] R. Aloisio and V. S. Berezinsky. Anti- 
GZK Effect in Ultra-High-Energy Cosmic 
Ray Diffusive Propagation. Ap. J., 625:249, 
2005. 

[7] R. Aloisio and V. S. Berezinsky. Diffu- 
sive Propagation of Ultra-High-Energy Cos- 



mic Rays and the Propagation Theorem. Ap. 
J., 612:900, 2006. 
[8] E. Amato and P. Blasi. A general so- 
lution to non-linear particle acceleration at 
non-relativistic shock waves. MNRAS Lett., 
364:76, 2005. 
[9] E. Amato and P. Blasi. Non-hnear particle 
acceleration at non-relativistic shock waves 
in the presence of self-generated turbulence. 
MA^/?A5, 371:1251,2006. 

[10] A. E. Bell. The acceleration of cosmic rays 
in shock fronts. I. MNRAS, 182:147, 1978. 

[11] A.E.Bell. The non-linear self-regulation of 
cosmic ray acceleration at shocks. MNRAS, 
225:615, 1987. 

[12] A. E. Bell. Turbulent amphfication of mag- 
netic field and diffusive shock acceleration of 
cosmic rays. MNRAS, 353:550, 2004. 

[13] E. G. Berezhko and D. C. Ellison. A Simple 
Model of Nonlinear Diffusive Shock Accel- 
eration. Ap. J., 526:385, 1999. 

[14] E. G. Berezhko and H. J. Volk. Spec- 
trum of cosmic rays, produced in super- 
nova remnants. In ICRC, Merida, Mexico, 
ICRC.OG.Olll, 2007. 

[15] E. G. Berezhko and H. J. Volk. Theory of 
cosmic ray production in the supernova rem- 
nant RX J 17 13. 7-3946. hi ICRC, Merida, 
Mexico, ICRC.OG.0614, 2007. 

[16] Puehlhofer G. Berezhko, E. G. and H. J. 
Volk. Cosmic ray acceleration and gamma 
ray production in the supernova renmant RX 
J0852.0-4622. In ICRC, Merida, Mexico, 
ICRC.OG.0597, 2007. 

[17] Gazizov A. Berezinsky, V. and S. Grigorieva. 
Dip in UHECR spectrum as signature of pro- 
ton interaction with CMB. Phys. Lett. B, 
612:147, 2005. 

[18] Amato E. Blasi, R and D. CaprioH. The 
maximum momentum of particles acceler- 
ated at cosmic ray modified shocks . MNRAS, 
375:1471,2007. 

[19] P. Blasi and E. Amato. A kinetic approach 
to non resonant modes and growth rates 
of streaming instability: consequences for 
shock acceleration. In ICRC, Merida, Mex- 
ico, ICRC.OG.0342, 2007. 

[20] P. Blasi and E. Amato. Theory of non lin- 
ear particle acceleration at shocks and self- 



30th International Cosmic Ray Conference 



generation of the magnetic field. In ICRC, 
Merida, Mexico, ICRC.OG.0341, 2007. 
[21] P. J. Boyle et al. Cosmic Ray Energy Spec- 
tra of Primary Nuclei from Oxygen to Iron: 
Results from the TRACER 2003 LDB FHght. 
In ICRC, Merida, Mexico, ICRC.OG.1192, 
2007. 

[22J R. Biiler et al. (HESS Coll.). Energy 
spectrum of cosmic iron nuclei measured 
by HESS. In ICRC, Merida, Mexico, 
ICRC.OG.0284, 2007. 

[23] S. Coutu et al. Design and performance in the 
first flight of the transition radiation detec- 
tor and charge detector of the CREAM bal- 
loon instrument. Nucl. Instr. and Meth. A, 
572:485, 2007. 

[24] C. De Donato and G. Medina Tanco. The end 
of the galactic spectrum. In ICRC, Merida, 
Mexico, ICRCOG.1249, 2007. 

[25] Blasi P. De Marco, D. and T. Stanev. Nu- 
merical Propagation of Cosmic Rays in 
the Galaxy. In ICRC, Merida, Mexico, 
ICRC.OG.0736, 2007. 

[26] Blasi P. De Marco, D. and T. Stanev. Numeri- 
cal propagation of high energy cosmic rays in 
the Galaxy: I. Technical issues. JCAP, 6:27, 
2007. 

[27] Mastichiadis A. Dimitrakoudis, S. and 
A. Geranios. Simulation of Cosmic Ray 
propagation in the Galactic Centre Ridge 
in Accordance with Observed VHE gamma 
ray Emission. In ICRC, Merida, Mexico, 
ICRC.OG.l 127, 2007. 

[28] Grasso D. Springel V. Dolag, K. and 
I. Tkachev. Constrained simulations of the 
magnetic field in the local Universe and the 
propagation of ultrahigh energy cosmic rays. 
JCAP, 1:9, 2005. 

[29] Van Eldik C. Hinton J. (HESS Collaboration) 
Egberts, K. Measurement of cosmic ray elec- 
trons with HESS. In ICRC, Merida, Mexico, 
ICRC.OG.0446, 2007. 

[30] J. J. Engelmann et al. (HEAO Coll.). Astwn. 
andAstrop., 233:96, 1990. 

[31] S. Funk. VHE Gamma-ray supernova rem- 
nants. astro-ph/0701471, 2007. 

[32] T. Hams et al. Results from BESS-Polar I 
2004 Antarctica Flight. In ICRC, Merida, 
Mexico, ICRC.0G.1119, 2007. 



[33] M. HiUas. Can diffusive shock acceleration in 
supernova renmants account for high-energy 
galactic cosmic rays? J. ofPhys. G, 31:95, 
2006. 

[34] C.-Y Huang and M. Pohl. Monte Carlo 
Study of Cosmic-Ray Propagation in the 
Galaxy and Diffuse Gamma-Ray Production. 
In ICRC, Merida, Mexico, ICRC.OG.1087, 
2007. 

[35] Waelkens A. H. Farrar G. R. Jansson, R. and 
T. A. EnssUn. Large scale magnetic field of 
tiie MilkyWay from WMAP3 data. In ICRC, 
Merida, Mexico, ICRC.OG.1271, 2007. 

[36] J. R. Jokipii. Cosmic-Ray Propagation. I. 
Charged Particles in a Random Magnetic 
Field. Ap. J., 146:480, 1966. 

[37] J. R. Jokipii and J. Giacalone. Effects of 
Large-Scale Upstream Turbulence on a Su- 
pernova Blast Wave. In ICRC, Merida, Mex- 
ico, ICRC.OG.0078, 2007. 

[38] Das S. Ryu D. Kang, H. and J. Cho. Propa- 
gation of UHE protons through a magnetized 
large scale structure. In ICRC, Merida, Mex- 
ico, ICRC.OG.04I8, 2007. 

[39] P. O. Lagage and C. J. Cesarsky. Cosmic- 
ray shock acceleration in the presence of self- 
excited waves. Astron. andAstrop., 118:223, 
1983. 

[40] R O. Lagage and C. J. Cesarsky. The 
maximum energy of cosmic rays accelerated 
by supernova shocks. Astron. and Astrop., 
125:249, 1983. 

[41] M. Lemoine. Extragalactic magnetic fields 
and the second knee in the cosmic-ray spec- 
trum. Phys. Rev. D, 71:3007, 2005. 

[42] Pelletier G. Lemoine, M. and B. Revenu. On 
the Efficiency of Fermi Acceleration at Rela- 
tivistic Shocks. Ap. J. Lett., 645:129, 2006. 

[43] M. A. Malkov. Analytic Solution for Nonlin- 
ear Shock Acceleration in the Bohm Limit . 
Ap. J., 485:638, 1997. 

[44] M. A. Malkov and L. O'C. Drury. Non Hn- 
ear theory of diffusive acceleration of parti- 
cles by shock waves. Rep. Prog, in Phys., 
64:429, 2001. 

[45] J. F. McKenzie and H. J. Volk. Non-hnear 
theory of cosmic ray shocks including self- 
generated Alfven waves. Astron. andAstrop., 
116:191, 1982. 



Origin and Propagation of Cosmic Rays 



[46] I. V. Moskalenko and A. W. Strong. Isotopic 
composition of cosmic ray sources. In ICRC, 
Merida, Mexico, ICRC.OG.0738, 2007. 

[47] J. Niemiec and M. Pohl. Magnetic turbu- 
lence production by streaming cosmic rays 
upstream of SNR shocks. In ICRC, Merida, 
Mexico, ICRC.OG.1047, 2007. 

[48J Ostrowski M. Niemiec, J. and M. Pohl. The 
inefficiency of the first order Fermi process 
in UHECR production at relativistic shocks. 
In ICRC, Merida, Mexico, ICRC.OG.I05I, 
2007. 

[49] A. D. Panov et al. Relative abundances of 
cosmic ray nuclei B-C-N-O in the energy 
region from 10 GeV/n to 300 GeV/n. Re- 
sults from ATIC-2 (the science ight of ATIC). 
In ICRC, Merida, Mexico, ICRC.OG.00I8, 
2007. 

[50] I. H. Park et al. SiUcon charge detector for the 
CREAM experiment. Nucl. Instr. and Meth. 
A, 570:286, 2007. 

[51] N. H. Park et al. Relative abundances of 
heavy ions measured by the CREAM-II sil- 
icon charge detector. In ICRC, Merida, Mex- 
ico, ICRC.OG.0593, 2007. 

[52] P. Picozza et al. The Physics of PAMELA 
Space Mission. In ICRC, Merida, Mexico, 
ICRC.OG.0289, 2007. 

[53] V. S. Ptuskin and V. N. ZirakashviU. On 
the spectrum of high-energy cosmic rays pro- 
duced by supemova remnants in the pres- 
ence of strong cosmic-ray streaming instabil- 
ity and wave dissipation. Astron. andAstrop., 
429:755, 2005. 

[54] B. F. Rauch et al. Measurement of the Rela- 
tive Abundances of the Ultra-Heavy Galactic 
Cosmic Rays (30 < Z < 40) with TIGER. 
In ICRC, Merida, Mexico, ICRC.OG.0I87, 
2007. 

[55] E. S. Seo et al. Cosmic Ray Energetics And 
Mass (CREAM) Overview. In ICRC, Merida, 
Mexico, ICRC.OG.0677, 2007. 

[56] A. Shalchi and R. SchUckeiser. Altemative 
explanation of the abundance ratio of sec- 
ondary to primary galactic cosmic ray nuclei. 
In ICRC, Merida, Mexico, ICRC.OG.0046, 
2007. 

[57] Miniati F. Sigl, G. and T. Ensshn. Ultra- 
high energy cosmic ray probes of large scale 



structure and magnetic fields. Phys. Rev. D, 
70:3007, 2004. 

[58] Kieda D. B. Swordy, S. P and S. P Wakely. A 
high resolution method for measuring cosmic 
ray composition beyond 10 TeV. Astropart. 
Phys., 15:287,2001. 

[59] J. Vink. X-ray high resolution and imaging 
spectroscopy of supemova remnants . In The 
X-ray Universe 2005, page 319, 2006. 

[60] Berezhko E. G. Volk, H. J. and L. T. Kseno- 
fontov. Magnetic field amplification in Tycho 
and other shell-type supernova remnants. As- 
tron. andAstrop., 433:229, 2005. 

[61] S. A. Wissel et al. The Status of the Track 
Imaging Cerenkov Experiment. In ICRC, 
Merida, Mexico, ICRC.OG.073I, 2007. 

[62] S. A. Wissel et al. (Veritas Coll.). Studies of 
Direct Cherenkov Emission with VERITAS. 
In ICRC, Merida, Mexico, ICRC.OG.0732, 
2007. 

[63] Y. S. Yoon et al. H and He spectra from the 
2004/05 CREAM-I flight. In ICRC, Merida, 
Mexico, ICRC.OG.0778, 2007. 

[64] H. S. Yoshida et al. Energy spectrum and 
arrival directions of high-energy electrons 
over 100 GeV observed with PPB-BETS. 
In ICRC, Merida, Mexico, ICRC.OG.0892, 
2007. 

[65] R. Zei et al. Preliminary measurements of 
Carbon and Oxygen energy spectra from the 
second flight of CREAM. In ICRC, Merida, 
Mexico, ICRC.OG.030I, 2007. 

[66] Pochepkin D. N. Ptuskin V. S. Zirakashvili, 
V. N. and S. I. Rogovaya. Propagation 
of ultra-high-energy cosmic rays in Galactic 
magnetic fields. Astron. Lett., 24:139, 1998. 



